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Abstract. The propagation of a weak probe field in a laser-driven four-level atomic 
system is investigated. We choose mercury as our model system, where the probe 
transition is in the ultraviolet region. A high-resolution peak appears in the optical 
spectra due to the presence of interacting dark resonances. We show that this narrow 
peak leads to superluminal light propagation with strong absorption, and thus by itself 
is only of limited interest. But if in addition a weak incoherent pump field is applied to 
the probe transition, then the peak structure can be changed such that both sub- and 
superluminal light propagation or a negative group velocity can be achieved without 
absorption, controlled by the incoherent pumping strength. 
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1. Introduction 

Optical properties of an atomic medium can be substantially modified by the application 
of external fields. In particular, atomic coherence induced by laser fields plays an 
important role in light-matter interaction and has found numerous implementations 
in optical physics [1]. One prominent application is the modification of the propagation 
of a light pulse through an atomic medium, which depends on the dispersive properties 
of the medium. The study of such pulse propagation phenomena has been triggered 
by a series of papers by Sommerfeld and Brillouin [2l [3] and continues to be of much 
interest [H El El [71 El [9] . It is well known that the group velocity of a light pulse can be 
slowed down [lOl [11], can become faster than its value c in vacuum, or can even become 
negative [121 [13] • Note that superluminal light propagation with group velocity larger 
than c cannot transmit information faster than the vacuum speed of light [H] , such that 
it is not at odds with causality. Superluminal light propagation has been investigated for 
many potential uses, not only as a tool for studying a very peculiar state of matter, but 
also for developing quantum computers, high speed optical switches and communication 
systems [TSj . 

Both experimental and theoretical studies have been performed to realize super- 
and subluminal light propagation in a single system. For example, speed control in 
atomic systems has been achieved by changing the frequencies, amplitudes or phase 
differences of the applied fields. It has been shown that switching from subluminal 
to superluminal pulse propagation can be achieved by the intensity of the coupling 
fields [IHl [T71 [HI [19] , and the relative phase between two weak probe fields [20]. Morigi 
et al. [21] have compared the phase-dependent properties of the o (diamond) four level 
system with those of the double A system. In Ref. [13], gain-assisted superluminal 
light propagation was observed in a cesium vapor cell while in most other studies, 
superluminal light propagation is accompanied by considerable absorption. Sub- and 
superluminal light propagation together with nonlinear optical gain or losses were 
observed in [22]. Two of the present authors suggested to use an incoherent pump 
field to control light propagation from subluminal to superluminal [23l [23] . Recently, 
we have studied the light propagation of a probe pulse in a four-level double lambda 
system, where the applied laser fields form a closed interaction loop [25] . In such systems, 
the finite frequency width of a probe pulse requires a time dependent treatment of the 
light propagation. We have found both sub- and superluminal light propagation without 
absorption or with gain, controlled by the Rabi frequency of one of the coupling fields. 

All these effects depend on the modification of the dispersive and absorptive 
properties of the atomic medium. A particular class of systems that allows to 
modify the optical response to a great extend are those with so-called interacting dark 
resonances [26] . A characteristic feature of such systems is the appearance of very sharp, 
high-contrast structures in the optical spectra. Resonances associated with double dark 
states can be made absorptive or transparent and their optical properties such as width 
and position can be manipulated by applying suitable coherent interactions. It was 
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Figure 1. (a) Energy scheme of the four level atomic system considered. Transition 
|2) ^ |4) is driven by a strong laser field, transition |1) <-> |4) by a weak coupling field, 
and the probe field interacts with transition |2) <-^- |3). In addition, a weak incoherent 
field is applied to the probe field transition, (b) A possible realization of the scheme in 
mercury. Population transfer to state (t'Pq has to be compensated via a repump field. 



also shown that very weak incoherent excitation of the atoms can be sufficient to turn 
absorptive features into optical gain structures. This has been proposed as a model 
system to obtain strong laser gain in the ultraviolet and vacuum ultraviolet regime by 
Fry et al. [27]. 

In this paper we consider probe pulse propagation through a system which exhibits 
interacting dark resonances. The level configuration of our four-level scheme is based 
on the lasing system proposed in |23, and consists of three atomic states in ladder 
configuration, with an additional fourth perturbing state coupled by a laser field to the 
upper state of the ladder system. The lower transition of the ladder system acts as 
the probe transition. This system can be realized, e.g., in mercury, where the probe 
transition has a low wavelength of 253.7 nm, i.e., in the ultraviolet region. We find 
that the medium susceptibility in dependence on the probe field detuning exhibits high- 
contrast structures characteristic of interacting dark states. These structures typically 
lead to superluminal probe field propagation with high absorption, and thus as such 
are of limited interest. If, however, a weak incoherent pumping is applied in addition 
to the probe field transition, then we find that in the region around a narrow structure 
both sub- and superluminal propagation as well as negative group velocities are possible 
without absorption, controlled by the incoherent pumping strength. 

2. Analytical considerations 

2.1. The Model System 

We consider an atomic four level system as shown in figure [T](a). Transition |2) ^ |4) 
is driven by a strong coherent field with frequency 0^42 and Rabi frequency (742. A 
weak coupling field with frequency u;4i and Rabi frequency (741 is applied to transition 
|1) ^ |4). The weak probe field with frequency UJ22, and Rabi frequency (723 = Qp couples 
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to transition |2) ^ |3). Finally, an incoherent driving field with pump strength A is 
applied to the probe transition. We further include spontaneous decay with rates 741, 
742, 723, and 7i3, respectively, on the dipole- allowed transitions. The atomic transition 
frequencies are denoted by Uij, and the laser field detunings with respect to the atomic 
transition frequencies are Ajj = ojij — COij {i,j G 1, ... ,4). A realization of our level 
scheme can be found, e.g., in mercury, see figure [T](b). 

The density matrix equations of motion, in the rotating wave approximation, are 

Pii = - 2713P11 + 2741P44 - igl^pu + ignPAi , (la) 
P22 = - 2723P22 + 2742P44 - 2Ap22 + 2Ap33 

+ i9*pP32 - igpP23 - i9l2p2i + ^fi'42P42 , (1^) 

P33 = 2713P11 + 2723P22 + 2Ap22 - 2Ap33 - iglpz2 + igpP2?. , (Ic) 
P12 = - (ri2 + iA4i - zA42 + A)pi2 - igl2Pi4 " WpPi3 + ^5'4iP42 , (Id) 
Pi3 = - (ri3 + iA4i - iA42 - iAp + A)pi3 - iglpi2 + ig^iP^is , (le) 

Pl4 = - (Fm + zA4i)pi4 - ig4lPn + «fi'4lP44 - ^fi'42Pl2 , (If) 

P23 = - 0^23 - i\ + 2A)p23 - ig*p22 + ig*pp33 + «fi'42P43 , (1^) 

P24 = - (r24 + ^A42 + A)p24 - ifi'42P22 + ^fl'42P44 + «5'pP34 " ^5'4lP21 , (1^^) 

P34 = - (r34 + i\ + ^A42 + A)p34 + igpP2A " «fi'4lP31 " «5'42P32 , (1«) 

P44 = 1 - Pll - P22 - P33 • (Ij) 

In the above equations, Fjj = (27j + 27j)/2 are the damping rates of the coherences with 
7j being the total decay rate out of state and Ap = A23 is the probe field detuning. 

Our main observable is the response of the atomic medium to the probe field. 
As will be discussed in Sec. 12. 2[ the linear susceptibility of the weak probe field is 
determined by the probe transition coherence P23. We therefore proceed by solving the 
above equations (|Taj)-(nij) in the steady state under the assumption of specific parameter 
relations. 

First, in the absence of the incoherent pump field (A = 0), an expansion of the 
steady state coherence P23 to the leading order in the probe field Rabi frequency gp 
yields 



-5'p(lfl'41 


P ^ Cl3 ■ ^^34) 


lfl'41 


^^^23 + C'l3( ^742 


P ^ C23 ■ C*34) 



^ 1(741 pC'23 + C'l3(|fi'42P ^ C23 ■ C*34) ' 

Ci3 = Ap - A41 + A42 + ^ Fi3 , (26) 

C34 = Ap + A42+^F34, (2c) 

C23 = Ap + «F23. {2d) 

It will turn out that an interesting parameter range for the present study is given by 

A41 = A42 = 0, (3a) 

Ap < 731,741,742 , (36) 

6-41 < gi2 , (3 c) 

713 =0. (3d) 
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In this limit, equation becomes 

P23 = zMi^ ^^El^ . (4) 



-gp{\gAi? - 


i\T^4) 


\gi2?^p + ^[|5'4iPr23 


- A2(r34 + r23)] 



An inspection of equation (jlj) reveals that the imaginary part is strictly positive, and 
the half width of the absorption peak around Ap = is determined by 

~ ( — ) 123 = ( — ) 723 • (5) 

Next, we seek the corresponding steady state solution for p23 with incoherent pump 
field with pump intensity A. The parameters are chosen to satisfy equations (!3aj) -( l3dll 
as well as the new condition on the pump field 

Ao < A < 741,742. (6) 

Further, we assume the Rabi frequencies gij to be real in the following. We obtain in 
leading order of the probe field coupling Qp 

^ 9li9pl23 Ap-iA 

(^?|2 723 + 2Ar24 742)A2 + A2- 

Here the parameter Aq is defined by 



941 723 (741 +723) f 941 



^0 = — T ^ 723 • (8) 

5'42 741 + 723 r34 (741 + 723) V fl'42 / 

Since |5'4i/fi'42p723 can be made small, for a suitable combination of the Rabi frequencies 
g^i and (742 the condition A ^ Aq can be fulfilled even for incoherent pump strengths 
which are orders of magnitude smaller those required, e.g., to saturate the optical 
transition. 

We find that the imaginary part of equation ([7]) is negative if the condition A Aq 
is fulfilled. Thus Aq indicates the incoherent pumping rate at which the absorption peak 
turns into a gain structure, if the conditions in equations (l3a|) -([3dl) and ([6]) are fulfilled. 

2.2. Observables 

Our main observable is the response of the atomic medium to the probe field. The linear 

susceptibility of the weak probe field can be written as [28] 

2Nri„ , , 

XM = —^P23{^p) , (9) 

where N is the atom number density in the medium, r]p is the probe transition dipole 
moment and X = x' + ^x"- The real and imaginary parts of xi^p) correspond to the 
dispersion and the absorption, respectively. The slope of the dispersion with respect 
to the probe detuning has a major role in the calculation of the group velocity. We 
introduce the group index, rig = c/vg, where the group velocity Vg of the probe field is 
given by [ini Hn] 

c 



1 + 27rx'M + 27rup^^ 



(10) 
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Equation (ITU!) implies that, for a negligible real part x'iy^p)) the group velocity can be 
significantly reduced via a steep positive dispersion. Strong negative dispersion, on the 
other hand, can lead to an increase in the group velocity and even to a negative group 
velocity. 

Substituting equations (JTj) and ([9]) in equation (ITOl) . the group index of the probe 
field evaluates to 

n -1 = ^'i^^^^^ (11) 

It can be expected from equation ffTTj) that for suitable parameters, the group index 
around Ap = is negative and accompanied by gain, and this is indeed what we find 
below. 

The relation between coherence and susceptibility equation ([9]) can be rewritten as 

where we have used 723 = {ripU}^^) / (2,TieQhc^) and = ri23Ep/h as well as UJ23 = 2itc/\23 
with the probe transition wavelength A23. For mercury probe wavelength 253.7 nm, 
particle density = lO^^cm"^ and 723/7 = 0.14 as found in mercury one finally obtains 

XK) = 1.74 X 10- ^ . (13) 

Throughout our discussion of numerical results, we will assume these parameters in 
order to evaluate the susceptibility. 

2.3. Dressed-state analysis 

We now introduce the dressed states generated by the strong driving field acting on 
transition |2) ^ |4) and the coupling field acting on transition |1) ^ |4), in order to 
demonstrate the presence of interacting dark resonances due to the perturbing field with 
Rabi frequency g^n . In the absence of the incoherent pump field, the dressed states 
are 

10) = - , |1) + ,f , |2) , (14«) 

^^{9li+9l2) ^^{9li + 9I2) 

with energies 



Ao = 0, \± = ±h.^gl^+gl2. (15) 

The two dressed states |±) correspond, in the limit of vanishing driving field (741, to 
the usual Autler-Townes dressed components split by 2%42. The third dressed state |0) 
coincides in this limit with the bare state |1) and hence is decoupled from the fields. This 
is no longer so in the presence of a second weak driving field g^. In this case the dressed 
state |0) contains an admixture of |2) and thus has a nonzero dipole matrix element 
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Figure 2. Real (blue dashed) and imaginary (red solid) parts of the susceptibility 
X as a function of the probe detuning Ap for the parameters 741 — 7, 723 = 0.14 7, 
742 = 0.797, 713 = 0.01 7, 9p - 10-^7, 541 = 0, 542 - 47, A = 0, A42 - A41 = 0. 

with the state |3). As a result of this couphng, there are transitions between |0) and 
|3), corresponding to three photon resonances from |1) to |3) that exhibit interference 
effects [27]. 

3. Results 

In figure [2] we show the real (blue dashed) and imaginary (red solid) part of the probe 
field susceptibility x versus the probe detuning Ap, which correspond to the dispersive 
and absorptive properties of the medium, respectively. In this figure, the perturbing 
laser field is switched off, 5^41 = 0. The other parameters are 741 = 7,723 = 0.14 7, 
742 = 0.797, 713 = 0.01 7, Qp = 10-^7,(/42 = 47, A = 0, A42 = A41 = 0. Note that 
the ratios of the decay rates correspond to the case found in mercury, see figure [D^b) . 
We have added a weak decay rate 713, since otherwise in the steady state all population 
is trapped in |3). The driving field with Rabi frequency (742 leads to an Autler-Townes 
doublet with a dip in the absorption at zero detuning, i.e., partial electromagnetically 
induced transparency (EIT). The slope of the real part of the susceptibility in the region 
of reduced absorption is positive. We thus find that subluminal light propagation occurs 
around zero detuning with reduced absorption as it is common for EIT. If the state 
1 4) was long-lived, then the EIT leading to the partial transparency would be more 
pronounced such that the absorption would vanish at zero detuning. 

In figure [3], in addition we apply the weak perturbing field with Rabi frequency 
5'4i = 0.047, and assume negligible decay on transition |3) ^ |1), since a trapping in 
this state is now avoided by the additional laser field. The results are identical to figure [2] 
except for a narrow absorption spike at around zero detuning. The shape and width of 
the absorption spike are determined by equation fl2aj) and equation ([5]), respectively. In 
particular, the width is much less than the natural linewidth. Again, for a long-lived 
state |4), the transparency regions on each side of the absorption spike would become 
two points of EIT, i.e., a double dark state [26]. In terms of the light propagation, the 
slope of the real part of the susceptibility around zero detuning is negative such that 
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Figure 3. Real (blue dashed) and imaginary (red solid) parts of the susceptibility x 
as a function of the probe detuning for 713 = and g^i — 0.04 7. The other parameters 
are the same as in figure [2] (b) is a closeup on the central part of (a). 




Figure 4. Real (blue dashed) and imaginary (red solid) parts of the susceptibility x 
as a function of the probe detuning for A = 4 x 10~^ 7. The other parameters are the 
same as in figure [3l (b) is a closeup on the central part of (a). The purple vertical 
lines indicate the roots of the imaginary part. 



superluminal light propagation could be observed, albeit with high absorption. 

We now in addition apply a weak incoherent pumping field on the probe transition 
|2) — |3). Figure m shows the corresponding results. The incoherent pump field rate 
is chosen as A = 4 x 10~^7. In this case, the superluminal light propagation found 
in figure [3] at Ap = switches to subluminal propagation, and the absorption spike at 
zero detuning becomes a gain spike. The shape of the gain spike is determined by the 
imaginary part of equation ([7]) which is Lorentzian with halfwidth equal to A. This can 
be understood from the fact that the spike arises from a three-photon transition from 
|1) to |3). In the absence of the incoherent pump field, most of the population is in level 
1 3). Therefore, the probe field is absorbed. But in the presence of the weak incoherent 
pump field, the population is transferred to level |1). This optical pumping leads to the 
gain spike in the spectrum. 

Furthermore, it can be seen from Fig. |4] that at Ap ^ ±3.1 x 10~^7 (indicated 
by the purple vertical lines), the imaginary part of the susceptibility vanishes together 
with a negative slope of the real part. At these probe field detunings, the real part of 
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Figure 5. (a) Imaginary part x" of the susceptibility as a function of the pumping 
field strength A. The parameters are as in figure [3] with Ap = 0, and (742 = 157 (solid 
red), IO7 (short-dashed blue), 47 (long-dashed green), (b) Slope of the real part x' of 
the susceptibility at zero probe field detuning Ap = for parameters as in (a). 



the susceptibility itself is non-zero, and is negative (positive) for Ap ^ —3.1 x 10~^7 
(Ap —3.1 X 10~^7). In the following, we discuss the two cases of interest with resonant 
or non-resonant probe field separately. 

We start with the resonant case Ap = 0. In figure [S^a), we study the effect of 
the incoherent pumping strength A on the magnitude of the imaginary part of the 
susceptibility x at resonance Ap = 0. It can be seen that depending on the coupling 
field Rabi frequency 5(42, the transition from absorption to gain occurs at different values 
of the incoherent pumping. For the parameters of figure [3], which correspond to the long- 
dashed green curve in figure[5](a), the transition is at about A ^ 2 x 10~^ 7. This explains 
why gain could be observed for the parameters in figure HI On increasing the incoherent 
pumping further, the imaginary part approaches zero again. 

After the discussion of the absorption, we now turn to a discussion of our main 
observable, the group velocity. Since the real part of x itself vanishes at Ap = 0, 
the group velocity is determined by the slope of the real part of the susceptibility at 
Ap = 0, see equation (fTOl) . This quantity is shown in figure [5](b). It can be seen that, for 
no or small incoherent pumping, the system exhibits a negative slope, which leads to a 
superluminal or even negative group velocity. On increasing A, the slope can be adjusted 
to large positive values, where subluminal light can be expected. Thus in principle the 
system allows for a wide range of group velocities, controlled via the incoherent pump 
rate A. But from a comparison of figures [5] (a) and (b) it can be seen that typically 
negative slopes are accompanied by absorption, while positive slopes occur together 
with gain. Thus at Ap ?s 0, only a reduction of the group velocity is accessible in 
experiments without absorption. The different curves in figure [5] further show that the 
precise response of the system to the incoherent pumping can be controlled by varying 
the coupling field Rabi frequency 5^42- In particular, stronger coupling fields (742 may 
be favourable, since then the range of possible slopes is increased, as can be seen from 
figure [5](b). 

We now turn to a discussion of the non-resonant case, Ap 7^ 0, and focus on the 
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Figure 6. (a) Position Ao of the root of the imaginary part of the susceptibihty x as a 
function of the pumping field strength A. At this frequency, Ught passes unattenuated 
through the medium, (b) Slope of the real part x' of the susceptibility at the detuning 
Aq with vanishing absorption. The parameters are as in figure [3] but with (742 = IO7 
(short-dashed blue), 342 = 77 (solid purple), and 342 = 47 (long-dashed green). 



regions with vanishing absorption, such as Ap ^ ±3.1 x 10~^7 in figure HI It can be 
seen that around these probe field detunings, the imaginary part of the susceptibihty 
vanishes, such that the probe field passes unattenuated through the medium. At the 
same time, the real part of the susceptibility is non-zero, and has a negative slope. 
Therefore, at these frequencies, superluminal or negative group velocities are accessible 
without absorption. In order to study this result in more detail, in figure EJ^a) we show 
the probe field detuning Aq at which the imaginary part of the susceptibility vanishes 
as a function of the incoherent pumping rate A. It can be seen that for no or small 
incoherent pumping A, there is always absorption such that no Aq can be found. Once 
A is large enough for a root in the imaginary part of the susceptibility to occur, the 
position of the root first increases rapidly with A, and then saturates. The required 
value of A also depends on the strength of the coupling field 5^42 as can be seen from 
figure [6](a). 

The corresponding figure Mjo) depicts the slope of the real part of the susceptibility 
as a function of A. It can be seen that by varying the pump field strength A, both positive 
and negative slopes can be achieved at frequencies where the medium absorption is zero. 
After passing through a maximum positive slope, the slope drops to a minimum negative 
slope and then slowly increases again towards vanishing slope. For every value of the 
coupling field Rabi frequency (742, optimum values of A can be identified where the 
slope is steepest and either positive or negative. The maximum absolute values of the 
slope are of order 10~^ s~^, such that the third term 27rujpdx' /du in the denominator of 
equation ( fTOl) for our probe transition varies between approximately —10^ and -|-10^. 
Therefore, strongly sub- and superluminal propagation as well as a large range of 
negative group velocities occur without absorption in our sample, controlled by the 
magnitude of the incoherent pumping. It should be noted that only very weak incoherent 
pumping is required, as can be seen from the scaling of the x-axes in figures [61 

Throughout this section, the figures [2H6] have been obtained from a numerical 
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Figure 7. Real (a,c) and imaginary (b,d) parts of the susceptibility x as a function 
of the probe detuning. The analytical results are shown as solid red lines, whereas our 
numerical results are shown as dashed blue lines. The parameters in (a,b) are as in 
figure [HI and equation Q is shown as the analytical result. The parameters in (c,d) 
are as in figure HI with equation ([7]) as analytical result. 

solution of the full density matrix equations ( ITffljl -inijl. In the following, we verify 
our approximate analytical expressions, equations fl2fflj) -fl7j). by a comparison to the 
exact numerical calculations. The result is shown in figure U\ where the solid red 
curves correspond to the approximate analytical solutions, whereas the blue dashed 
curves represent our numerical results. The approximate result equation (jl]) for the 
case without incoherent pump field is shown in comparison to the numerical data 
in figure [71^a,b). It turns out that in this case, the results from equation (jl]) are 
virtually identical to the corresponding numerical results. Equation for the case 
with incoherent pumping is compared to the numerical results in figure [7](c,d). Here, 
the analytic results only describe the qualitative behavior of the curves. The reason 
for this is that in this figure, we chose parameters for which the condition A ^ Aq in 
equation ([6]) is not well satisfied. If the incoherent pumping A is increased, the agreement 
of the approximate results with the numerical calculation improves. Thus we conclude 
that our analytical results describe the system well enough to allow for an optimization 
of the parameters towards a desired peak structure, as long as the conditions on the 
parameters are satisfied. 
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4. Conclusion 

We have discussed the dispersive and absorptive properties of a four-level atomic medium 
that exhibits interacting dark-state resonances. In our numerical analysis, we have 
focused on mercury atoms with an ultraviolet probe field wavelength of 243.7 nm. 
Due to the interacting resonances, a high-resolution structure appears in both the 
absorption and the dispersion spectra. A weak probe field tuned to this resonance 
usually experiences superluminal propagation with absorption. But if in addition a 
weak incoherent pump field is applied to the probe transition, then the superluminal 
light propagation changes to subluminal light propagation accompanied by no absorption 
or gain. Slightly off resonance, the probe field experiences a vanishing imaginary part 
of the susceptibility. At these off-resonant frequencies, the real part of the susceptibility 
itself is non-zero and has a slope depending on the incoherent pumping strength. Thus 
both sub- and superluminal light propagation as well as negative group velocities can 
be achieved without absorption. The control via the incoherent pump fields suggests 
potential applications, e.g., in optical switching devices or in controllable pulse delay 
lines for the ultraviolet frequency region. 
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